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Abstract. When applied to the blast wave formed by the ex- 
plosion of a massive star as a supernova (SN), the theory of 
diffusive particle acceleration at shock fronts predicts a very 
high energy density in cosmic rays. Almost immediately af- 
ter particles begin to be injected into the process, the cosmic 
ray pressure rises until comparable to the ram-pressure encoun- 
tered by the shock front. Those supernovae which are observed 
in the radio band i.e., radio supernovae (RSNe), provide direct 
evidence of particle acceleration in the form of synchrotron 
emitting electrons. Furthermore, these objects are particularly 
interesting, since they are usually surrounded by a relatively 
dense confining medium. The acceleration of cosmic rays can 
then lead to the production of very high energy ( VHE) gamma- 
rays which arise from collisions between energetic particles and 
target nuclei. We estimate the cosmic ray energy density assum- 
1-H ing a fraction ^ < 1 of the energy available at the shock front is 
■^ converted into cosmic rays. Combining this with the parameters 
^^ describing the environment of the SN progenitor, as deduced 
(~| from observations, and from more detailed modelling, we com- 
Q^pute the flux at Earth F(> 1 TeV) of photons of energy above 
Q 1 TeV. For the relatively weak but nearby supernova SN1987A 
5—1 we predict F{> 1 TeV) = 2 x 10~'^ photons s~' cm~^ before 
c/2 the shock front encounters the ring of dense matter seen by the 
^ Hubble Space Telescope. Subsequently, the flux is expected to 
rise further. The medium around SN1993J in M81 is thought to 
have a density profile p oc r~^l^, (with r the distance from the 
point of explosion) for which we predict a roughly constant flux 
of F(> 1 TeV) = 2 X 10~'^ photons s~' cm~^. Once the shock 
emerges into the region where p oc r~^, which is expected at 
larger distances, the flux should decrease in proportion to t~ ' . 
This object thus presents an interesting target for observation 
by telescopes which detect the Cerenkov light emitted by the air 
showers from VHE photons. A detection would provide obser- 
vational confirmation of cosmic ray acceleration in supernovae. 
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1. Introduction 

The strong shock fronts associated with supernova remnants 
are currently thought to be responsible for the acceleration of 
galactic cosmic rays of energy below about 10''* eV Apart from 
indirect arguments based upon the total energy budget required 
to sustain these particles against loss out of the galaxy, the 
main support for this hypothesis comes from the existence of a 
fairly well-developed theory for particle acceleration - diffusive 
acceleration - which has been successfully applied to super- 
nova shocks by several groups (Drury et al. 1989; Dorfi 1990; 
Jones & Kang 1990; Berezhko et al. 1994). However, there has 
not yet been any direct empirical confirmation of the theory. 
Estimates of the flux of high energy photons and neutrinos 
from galactic type II supernovae have been made by Berezin- 
skii & Ptuskin (1988). More recently, it has been pointed out 
that supernova remnants might be observable as sources of TeV 
photons using the atmospheric Cerenkov technique (Drury et 
al. 1994). Such an observation would be an important test of the 
theory, but the identification of a suitable target is not straight- 
forward (Aharonian et al. 1994). In this Letter, we show that 
the theory of diffusive particle acceleration predicts TeV pho- 
tons not only from supernova remnants (SNRs), but also from 
powerful radio supernovae shortly after explosion. 

Very high energy gamma-rays are produced by energetic 
protons (cosmic rays) when they interact with target nuclei, and 
create neutral pions. Thus, in any given volume of space, the 
production rate is proportional to the density of cosmic rays and 
the total mass of target nuclei. Compared to freshly exploded 
stars, supernova remnants occupy a relatively large volume con- 
taining a high target mass (of the order of a few solar masses). 
However, the pressure (or energy density) of cosmic rays within 
SNRs is much smaller than can be expected in the early phases 
of a supernova explosion which is confined inertially by a very 
dense medium (e.g., the stellar wind of the progenitor). Conve- 
niently, confined explosions also manifest themselves as radio 
sources as soon as the shock front has penetrated sufficiently 
far into the surrounding plasma to enable the photons to es- 
cape through it. Over a dozen such 'Radio Supernovae' (RSNe) 
have been observed to date. Of these, two are in galaxies close 
enough to our own to permit, in principle, detection of the 
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gamma-rays: SN1987A, which is an intrinsically weak source, 
but very close-by, and SN1993J, a very powerful RSNe situated 
in M81 at a distance of about 3.6 Mpc. 

2. Diffusive particle acceleration in SN1987A 

Despite the fact that theoretical arguments (Berezinskii & 
Ptuskin 1988; Volk & Biermann 1988) have emphasised that 
the early phase of a supernova explosion in an inhomogeneous 
medium is important for the acceleration of high energy cosmic 
rays, almost all detailed investigations of diffusive particle ac- 
celeration in SNRs have dealt with the case of a star exploding 
into uniform homogeneous surroundings. Furthermore, atten- 
tion has been focussed on the nuclear component of the accel- 
erated particles, since it is these which must carry off most of 
the energy and thus ultimately determine the overall accelera- 
tion efficiency. The case of SN1987A changed the situation in 
both respects. The progenitor star is known to have had a strong 
wind, which significantly modified the environment, and the 
theory of diffusive acceleration was applied in an effort to ex- 
plain the relativistic electrons responsible for the radio emission, 
rather than the acceleration of cosmic rays (Ball & Kirk 1992, 
henceforth BK). Although it is possible to construct models 
which account for the radio emission, the observed spectrum 
requires a shock of compression ratio about 2.7, much less than 
the value 4 expected for a strong supernova shock front ex- 
panding into the stellar wind of its progenitor. In view of the 
speed of the shock front (sa 20, 000 km s~'), it is unrealistic 
to suppose the Mach number could be reduced to the required 
value (about 2) merely by preheating the upstream medium. 
Instead, BK suggested that the shock front is modified by cos- 
mic rays, which exert a substantial pressure on the upstream 
medium. This kind of modified shock front is well-known in 
the nonlinear theory of diffusive shock acceleration, and can 
result in an extended upstream precursor and a gas subshock of 
relatively low compression ratio (e.g., Drury 1983). Following 
up on this suggestion, 'two-fluid' simulations similar to those 
of supernova remnants (e.g., Dorfi 1990) were performed for a 
supernova exploding into an inhomogeneous region. The first 
density profile to be investigated was 



ratio of 4, and that the flow remains spherically symmetric, we 
have: 



p(r) = M/(4--Kr^v„) , 



(1) 



corresponding to that of a progenitor undergoing steady mass- 
loss at a rate M via a stellar wind of constant velocity v^ 
(Duffy et al. 1993; Kirk et al. 1994a). Once the hydrodynamic 
structure of the flow was found, the radio emission was com- 
puted by considering the acceleration of electrons, which can 
be treated as test particles because of their negligible contri- 
bution to the energy density. More recently, the computations 
have been extended to a more realistic picture of the surround- 
ings of SN1987A (Duffy et al. 1994 henceforth DBK), which 
is thought to consist of a freely expanding wind of the form 
given in Eq. (1) up to a distance ?■„ from the explosion site, 
surrounded by a 'stagnation zone' of roughly constant density. 
Assuming the transition between the freely expanding wind and 
the stagnation zone is a strong shock front with a compression 



p(r) = M/(Trrlv„) for r > r^ 



(2) 



The stagnation zone, in turn, is presumably confined within 
a much denser, slower wind which flowed from the progenitor 
during its red-giant phase of evolution, up to about 10, 000 years 
ago. It is possible that the ring of dense material seen at about 
7" = 6 X 10'^ cm may form part of this wind (see McCray 1993 
for a review). 

A characteristic of all these calculations is that the shock 
front is strongly modified soon after cosmic ray injection is 
assumed to begin. This can be understood by considering the 
relationship between the expansion timescale and the timescale 
for the rise in cosmic ray pressure in the test-particle approx- 
imation (DBK). Nonlinear effects quickly become important 
and modify the shock once the pressure in cosmic rays be- 
comes comparable to the ram-pressure of the medium entering 
the shock front. An equivalent viewpoint is that this saturation 
point is reached when the cosmic rays take a substantial fraction 
^ < 1 of the kinetic energy of the matter flowing into the shock 
front from upstream. Writing 7"s(t) for the radius and Vs(t) for 
the velocity of the modified shock front in the rest frame of 
the progenitor at time t, we can estimate the energy density in 
cosmic rays Ec as 

Ec(t) ~ ^/ dt'47r[7-,(t')]' ^pKs(i')] [v.(t')]' , (3) 

where V(t) is the volume occupied by the cosmic rays, roughly 
equal to that of the shocked wind. We define a filling factor / 
as follows 



V(t) = f4w[v,(t)f/3 . 



(4) 



If pcomp is the overall compression ratio of the modified shock 
consisting of both sub-shock and precursor, then 1 > / > 

1 - (1 - 1/peomp)^ 

An unsolved question in the theory of diffusive acceleration 
which is important for the determination of the efficiency (j) 
is that of how particles are injected into the mechanism. The 
parameterisation of this process used in the two-fluid models 
of SN1987A is approximately the same as that used for models 
of supernova remnants. It is also similar to the more detailed 
kinetic models of supernova remnants presented by Berezhko et 
al. (1994). However, independent of the assumptions concern- 
ing injection, all treatments find a cosmic ray energy density in 
agreement with Eq. (3) once the shock is significantly modified. 

The basic quantity required for a calculation of the gamma- 
ray production rate is the product of the target mass Mtarget and 
the energy density Ec of cosmic rays. The target mass is given 
by the mass in the shocked wind: 
t 

(5) 



Mtarget(<) 



= / dt'4w[v,(t')fp[v,(t')]v,(t') . 
Jo 



During the free-expansion phase of the supernova, it is a good 
approximation to treat v^ as constant, in which case 



ScMtarget = 3cj)Mi,,^JifSwvd^) 



(6) 
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Drury et al. (1994) have given a convenient method of estimat- 
ing the gamma-ray production rate Nph given the target mass 
and cosmic ray energy density: 



Nph = g-E'cMtarget/M 



(7) 



where /x is the mean molecular weight of the progenitor wind 
and the quantity q depends on the spectrum of the cosmic rays 
and the energy of gamma-rays considered. For a wind profile of 
the formp oc t""", we have Mtarget tx t^"" and NphO^ t^"^". In 
the vicinity of the shock, the cosmic ray spectrum is expected 
to be relatively hard (Berezhko et al. 1994), and we therefore 
take q = 10~'^ for gamma-rays of energy greater than 1 TeV 
(Drury et al. 1994). Assuming a distance of 55kpc, the flux 
from SN1987A whilst it was expanding into the region r < r^, 
where the wind corresponded to a constant mass-loss rate is 
given by 



F(> 1 TeV) 



2x10-'^^ 



M 



Vs 



2 X lO^kms" 
photons s~' cm~ 



/ \ 10-5M@yr-i 

2 



500 km s" 




(8) 



This is in rough agreement with the estimate for a galactic type 
II SN by Berezinskii & Ptuskin (1988), if we take into account 
the different wind speed and greater distance to SN1987A. 
According to the best-fit model of DBK, 7"w = 2 x 10'^ cm, 
and the shock entered the stagnation zone in 1990. This is in 
agreement with the increase in soft X-ray emission observed 
at about this epoch (Beuermann et al. 1994). It follows that at 
present the gamma-ray flux is given approximately by (Kirk et 
al. 1994b): 



F{> ITeV) = 1.7 X 10 



-13^ 



M 



V, 



2 X 104 kms 
2 X 10'^ cm 



/ \ 10-5M@yr- 

5 / X -2 



-4.^.3 



500 km s" 

photons s~' cm~^ 



8yr 



(9) 



This flux level is very sensitive to the shock speed, which we 
expect to decrease somewhat in the stagnation zone. It also 
depends on the assumption of spherical symmetry, which is 
certainly inappropriate at distances approaching that of the ring. 
Once the dense ring is encountered we may expect the target 
mass to increase rapidly. Assuming, conservatively, that the 
cosmic ray pressure remains at the level in the stagnation zone 
and that the entire 0.05 Mq of the ring is engulfed, we find 



13^ 



F(> ITeV) = 3.3 x 10"'"^ 



p(r > 7-w) 
10-22gcm-3 



2 X 104 km s" 



— I photons s ' cm 



(10) 



3. SN1993J and other RSNe 

The radio emission of several supernovae has been successfully 
modelled as synchrotron radiation from relativistic electrons 
which move out with the shock front through the medium sur- 
rounding the progenitor (Chevalier 1982; Weiler et al. 1986). 
As the layer of absorbing material between the observer and the 
shock front gets thinner, the radio emission increases in strength, 
emerging first at high frequencies, and later at the more strongly 
absorbed low frequencies. The rising phase of the emission thus 
yields information about the density of the progenitor wind, and 
the speed at which the shock moves through it. 

This approach does not include a theory of particle accel- 
eration per se; it is simply assumed that a fixed fraction of the 
energy available at the shock front is converted into both mag- 
netic field and relativistic electrons. The distribution function of 
the electrons is assumed to be a power-law. Detailed modelling 
of these sources using the diffusive acceleration mechanism is 
still in progress (Ball & Kirk 1995) but it is sufficient for the 
present to note that in every case the observations require a 
steeper power-law spectrum than is predicted by diffusive ac- 
celeration at a strong shock. This suggests that in every RSNe 
the shock front is modified by the pressure of cosmic rays, so 
that the energy density in these can be estimated from Eq. (3). 

SN1993J is a particularly interesting RSN which is only 
3.6 Mpc distant and appears to have been a late-type supergiant 
with a massive, slow stellar wind. According to Van Dyk et 
al. (1994), the time dependence of the free-free absorption of 
radio emission by the wind in front of the supernova shock 
indicates that the mass-loss rate M of the progenitor's wind was 
not constant, but rather decreased as the star evolved towards 
explosion. This conclusion is supported by the behaviour of 
the X-ray emission (Fransson et al. 1994). The density profile 
around the progenitor is believed to be p oc r~^l^, from which 
it follows that 



M, 



target 



{v) = 2r M (r)/(3vj 



(11) 



where M (r) refers to the mass-loss rate at a time r/v^, before 
explosion. Van Dyk et al. (1994) find M= 1.5 x lO-^Mgyr"' 
at 7" = 4. 1 X lO'*" cm, and thus the flux is 



F(> ITeV) = 1.8 X 10 



Vs 



2x 104 kms"' 
photons s~' cm~^ 



,2^ / M (4.1 X 10"'cm) 
7 yi.Sx 10-4M@yr-i 

2 / X -2 



lOkms" 



(12) 



independent of the time of observation. Of course, the p oc 
r~^/^ density profile cannot extend indefinitely, since this would 
imply an infinitely large mass-loss. This suggests that the t~' 
behaviour of the constant M case (Eq. 8) should apply once the 
shock reaches a sufficiently large radius. The transition between 
these phases should also be evident from the radio data. 
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4. Discussion 

From the point of view of particle acceleration theory, the most 
uncertain factors in the estimated gamma-ray fluxes are the 
efficiency cj) and the filling factor /. The latter depends on the 
total compression ratio of the subshock plus precursor system 
and the extent to which the post-shock flow is decelerated. A 
total compression of 4 and constant rate of expansion leads to 
/ = 0.58. However, the true value of / may be significantly 
smaller, since not only do the cosmic rays soften the equation 
of state of the gas, they also remove energy from behind the 
shock by diffusion and both of these effects tend to increase 
Pcomp- The two-fluid simulations of SN1987A by DBK suggest 
typical overall compression ratios of pcomp ~ 10, and have 
/ K 0.27. 

An estimate of the ratio (j)/ f can be obtained directly from 
simulations, since for an t""" density profile and constant shock 
speed it is simply related to the ratio of the cosmic ray pressure 
Pc to the ram pressure ahead of the shock front: (j)/ f = 2(3 — 
n)(Pc/pvj)- In the models presented by DBK, this quantity rises 
during the free wind phase, and stabilises at a value of roughly 
Pel pvl K, 0.5, corresponding to an efficiency (j) ~ 30%. It 
follows that the flux of gamma-rays of energy above 1 TeV from 
SN1987A, as predicted in Eq. (10), may become detectable, 
when the shock encounters the dense ring known to surround 
the supernova. Furthermore, Eq. (12) suggests that SN1993J 
may be detectable, since it is only a factor of five lower than 
the flux from the Crab Nebula - the best observed TeV source 
to date. 

Such a high gamma-ray flux has direct implications for 
radiation in other bands. Electrons produced by the decay 
of TT"^ mesons created in cosmic ray-nucleon interactions in 
SN1993J have a luminosity which can be estimated, to or- 
der of magnitude, as 50 % of the gamma-ray luminosity: 
ie ~ 2 X lO^^ergs"'. Such electrons cool predominantly by 
synchrotron radiation and produce luminosities in the X-ray 
bands of ROSAT, ASCA and OSSE which are typically a factor 
of a few below L^. This estimate is not too far below observed 
fluxes (Zimmermann et al. 1994; Tanaka et al. 1993; Leising 
et al. 1994). Clearly, an accurate calculation of the relationship 
between the yield of gamma-rays above 1 TeV and the expected 
X-ray flux would be very interesting. This, however, requires 
detailed knowledge of the cosmic ray spectrum, which is not 
provided by two-fluid simulations, but must await a kinetic 
treatment, such as that performed for the wind of SN1987A by 
Berezhko (1994). 

Although we expect the two supernovae discussed above to 
be somewhat weaker than the brightest galactic SNRs (Aha- 
ronian et al. 1994) they have the advantage of being point 
sources for Cerenkov telescopes, and may thus be easier to 
observe. The detection of very high energy gamma-rays from 
either a supernova or a supernova remnant would provide strong 
support for the hypothesis that supernova remnants are respon- 
sible for accelerating galactic cosmic rays below lO'^^eV, and 
would also be a valuable diagnostic of the acceleration mecha- 
nism. 
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